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The thermodynamic properties of a one-dimensional model describing spin dynamics in the pres- 
ence of a twofold orbital degeneracy are studied numerically using the transfer-matrix renormaliza- 
tion group (TMRG). The model contains an integrable S'(7(4)-symmetric point and a gapless phase 
which is SC/(4) invariant up to a rescaling of the velocities for spin and orbital degrees of freedom 
which allows detailed comparison of the numerical results with conformal field theory. We pay spe- 
cial attention to the correlation lengths which show an intriguing evolution with temperature. We 
find that the model shows an intrinsic tendency towards dimerization at finite temperature even if 
the ground state is not dimerized. 
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I. INTRODUCTION 

Contrary to a usual band insulator the internal de- 
grees of freedom remain still active in a Mott insula- 
tor. In most of these insulators the spatial distribu- 
tion of electrons around every atom is frozen in at the 
melting point and changes little down to zero temper- 
ature so that the spin is the only active internal de- 
gree of freedom. This leads to the widely studied spin 
models as for example the Heisenberg model or the t-J 
model when hole-doping away from the insulating case 
is studied. In recent years, however, there has been con- 
siderable interest in systems which have additional or- 
bital degrees of freedom so that the spatial distribution of 
electrons (orbital order) may change considerably with 
temperature^^ Transition-metal oxides where orbitals 
are believed to play an important role are for example 
LaMnO^, LaTiOj 5 . 6 i 7 i 8 i 9 or YVOj 10 ! 11 - 12 ! 13 . 14 ! 15 . In 
such systems a large number of states may be nearly de- 
generate leading to very unusual magnetic properties as 
for example colossal magnetoresistanc& or magnetization 
reversals with temperatur o 11 : 12 but making on the other 
hand a microscopic description extremely difficult. 

From a theoretical point of view most interest has fo- 
cussed on one-dimensional (ID) versions of the Kugel- 
Khomskii Hamiltonian^ which allow for a full quantum 
mechanical treatment. The simplest and most symmetric 
model describing an 5 = 1/2 spin system with a twofold 
orbital degeneracy is derived if anisotropies in the orbital 
sector which are present due to Hund's rule coupling are 
completely ignored. In this case the Hamiltonian is given 
by 

H = E (SiSi+i + (r t T l+1 + ~) (1) 

where r = 1/2 is an orbital pseudospin describing the oc- 
cupation of two degenerate orbitals. The model has not 
only the obvious SU(2) x SU(2) symmetry, i.e. SU{2) 
symmetry in both spin and pseudospin space, but also 
5/7(4) symmetry unifying the spin and orbital degrees of 
freedomiiSiii The model is equivalent to the integrable 



q — 4 Uimin- Sutherland model, which has been solved 
by Bethe ansatz^ The exact ground state and the exci- 
tation spectrum consisting of three gapless mixed spin- 
orbital modes has been obtained. Thermodynamic prop- 
erties of the 5[/(4) symmetric model in ID have been 
investigated by the quantum Monte Carlo method^ In 
a real compound, however, symmetry will always be re- 
duced due to Hund's rule coupling. A widely studied toy 
model to investigate the physical consequences of sym- 
metry breaking is given by 

H = J2(S t S t+ i+x){T l T l+1 +y) (2) 

i 

with real constants x, yi 16 i 17 i 19 i 20 i 21 i 22 i 23 The model is 
57(2) x 57(2) symmetric and exhibits an additional Zi 
symmetry, interchanging spin and orbital degrees of free- 
dom, if x = y. The ground-state phase diagram depend- 
ing on x, y shows 5 phasesi£^iiSSiSi: a phase with fully 
polarized ferromagnetic spins and orbitals (I), phases 
with antiferromagnetic spin and ferromagnetic orbital 
states (II) and vice versa (III), a dimerized phase (IV) 
and a gapless phase (V) which shows 57(4) symmetry 
at low energies up to a rescaling of the spin and orbital 
velocities (see Fig. QJ. 

In this paper we will study thermodynamic properties 
of the Hamiltonian at the 57(4)-symmetric point, 
at x = 0.1, y = 0.25 and at x = 0.5, y — 0.5 which 
are representative for the gapless phase (V) and for the 
dimerized phase (IV), respectively. We will focus on two 
aspects: First we will show that all thermodynamic data 
at low temperatures (low-T) are consistently described 
by conformal field theory (CFT) if the system is criti- 
cal at T = 0. Second we will show that the model is 
highly susceptible for dimerization at finite temperature 
even if the ground state is not dimerized. The outline 
is as follows: In section [H] we will briefly introduce the 
transfer-matrix renormalization group (TMRG) which is 
the numerical method used to study the model at finite 
temperatures. We will discuss our results for the 57 (4)- 
symmetric case in section ITTT1 before we turn to the point 
x = 0.1, y = 0.25 (gapless phase) in IIVI and to x = 0.5, 
y = 0.5 (dimerized phase) in [V] A short summary and 
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FIG. 1: Phase diagram for the model ||2J taken from Ref. ITU 
Five phases exist as described in the main text. 



some conclusions are given in section fvTI 



II. TRANSFER-MATRIX RENORMALIZATION 
GROUP 

The idea of the TMRG is to express the partition func- 
tion Z of a one-dimensional quantum model by that of an 
equivalent two-dimensional classical model which can be 
derived by the Trotter-Suzuki formulai^^ For the classi- 
cal model a suitable transfer-matrix can be defined which 
allows for the calculation of all thermodynamic quanti- 
ties in the thermodynamic limit by considering solely the 
largest eigenvalue of this transfer matrix. Here we will 
use the Trotter- Suzuki mapping introduced in R.efs. l26l27l 
which yields 

Z = lim Tr{[Tx(e)T 2 (e)] M/2 } 
with 7i, 2 (e) = T R , L exp [-eH + 0(e 2 )] (3) 

where e = @/Mj& /? being the inverse temperature and 
M an integer number. Tr.£ denotes the right- and left- 
shift operator, respectively. This mapping allows for the 
definition of a transfer matrix with repeat length 1 and 
the free energy in the thermodynamic limit is then given 
by 



M 



-Tin A 



o 



(4) 



where An is the largest eigenvalue of the column-to- 
column transfer matrix Tm- This eigenvalue is non- 
degenerate and real for any finite temperature. Corre- 
lation lengths (CLs) £ n are also easy to calculate within 
this approach. They are given by 



C 1 = In 



Ao 
A„ 



= arg ' A^ 



(•5) 



where A„ are next-leading eigenvalues of Tm and k n is the 
corresponding wave vector. The appearance of in prin- 
ciple infinitely many CLs can be understood as follows: 
Every two-point correlation function can be expanded as 



(Oi T ) - (Ox) (O r ) = ]T M„ e-^» e h 



(6) 



with matrix elements M n where Oi denotes some opera- 
tor at site i. The long distance behavior is dominated by 
the CL £ Q belonging to the largest eigenvalue A Q (a ^ 0) 
which satisfies the condition M a ^ 0. Note that several 
CLs £ with the same wave vector k can appear in this 
asymptotic expansion. In the structure factor each expo- 
nential term yields an (in principle) measurable Lorentz 
function with center at k n , height ~ M n S, n /n and width 
~ 2/£„. The sharpest peak corresponds to the leading in- 
stability towards the onset of long range order and hence, 
a crossover in the leading CL indicates a change of the 
nature of the long range order. The transfer matrix is 
enlarged in imaginary time direction (being equivalent 
to a decrease in temperature) by using a density matrix 
renormalization group (DMRG) algorithm. In all fol- 
lowing calculations we have always retained between 350 
and 512 states. For details of the algorithm the reader is 
referred to Refs. 



III. THE SU(4) SYMMETRIC POINT 

Although thermodynamic properties of the SU (4) 
symmetric model have already been studied using quan- 
tum Monte Carlo algorithms^ we will start with this 
special point for several reasons. First, it provides a 
good test for the numerics because a comparison with 
Bethe ansatz and CFT results is directly possible. Sec- 
ond, there are still interesting properties of this model at 
finite temperature which have not been addressed so far: 
By comparing the numerically calculated specific heat, 
spin susceptibility and the correlation length data we will 
show that even in this highly symmetric case the system 
is susceptible for dimerization at finite temperature. We 
will also show that CFT describes consistently all data 
in the low temperature limit. 

In Fig. |2 the free energy / calculated by the TMRG 
method is shown. The low-energy properties of the model 
are known to be described by the universality class of the 
£17(4)1 Wess-Zumino-Witten (WZW) mo del equiv- 
alent to 3 free boson modes so that we expect from CFT 
in the low-T limit 



f = e 0- T T 
bv 



(7) 



with central charge c — 3. The velocity of the elemen- 
tary excitations v — tt/S and the ground state energy 
e w -0.20628 are known from BAi^ Eq. (J7J is in 
excellent agreement with the numerics for T/J < 0.1. 
The inset of Fig. [21 displays the entropy S = —df/dT 
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FIG. 2: Free energy (solid line) and low-T asymptotics 
(dashed line) / ~ eo — 4T 2 as expected from CFT with eo 
as specified in the text. The inset shows the entropy (solid 
line) with dashed lines denoting the low (S ~ 8T) and high-T 
asymptotics (S ~ In 4). 



which should be given in the low-T limit by S ~ ST 
according to Eq. JJJ and by S ~ In 4 at high temper- 
atures also in good agreement with the numerics. The 
specific heat calculated by taking the numerical deriva- 
tive C = -Td 2 f/dT 2 is shown in Fig. El Somewhat 




FIG. 3: Specific heat as a function of temperature (solid line) 
where the dashed line denotes again the low-T asymptotic 
behavior C ~ 8T. In the inset the spin susceptibility is shown 
and the dot denotes the value expected from CFT. 



within TMRG. 28 To calculate Xs(T = 0) within CFT one 
has to remind that the 577(4) i WZW model is equiva- 
lent to a sum of 2 decoupled SU(2) 2 WZW modelsiLSi 
This model is therefore a rare example of a critical spin 
model with Kac-Moody central charge k > 1. It has been 
conjectured^ that the susceptibility of an isotropic spin- 
s antiferromagnet will be given by Xs — k/2m) s so that 
we expect Xs = I/ttKs = 8/ir 2 which is in good agree- 
ment with a linear extrapolation of the numerical result 
to T = 0. 

By comparing C and \s which show both a peak in 
the same temperature region T/J ss 0.08 we could iden- 
tify this peak as due to the elementary excitations of the 
SU(A)x WZW model. Note also that this peak marks 
the point in temperature where the asymptotic expres- 
sions from CFT shown in Figs. |21 and |3| become invalid. 
To understand the appearance of the shoulder we note 
that the energy for a single bond is given by — J/4 if the 
spins form a singlet (triplet) and the orbitals a triplet 
(singlet) and by + J/4 if both are in singlet or triplet 
configuration. Therefore the excitation energy 8 for a 
single bond is 8 = J/2. In thermodynamic data, how- 
ever, it makes no difference for unbound excitations if 
they appear in pairs or not so that the energy scale is 
set by 8/2 = J/4. This approximately coincides with the 
shoulder in Fig. [3] So we may attribute this structure 
to the excitations of a chain with short range dimer or- 
der in both spin and orbital sector which seems to form 
at finite temperature although the ground state is not 
dimerized. This picture is further supported by an inves- 
tigation of the leading CLs which are plotted as T£/J in 
Fig. 0] because a 1/T divergence at low-T of those CLs 
belonging to the critical modes is expected (see below). 
Here a CL £d is present (dot-dashed line) which appears 
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surprisingly we find apart from a pronounced peak at 
T/J fa 0.08 an additional shoulder at T/J « 0.22. We 
will return to this point soon but first discuss the spin 
susceptibility x^ a t magnetic field h = which is shown 
in the inset of Fig. |3| This quantity is calculated nu- 
merically by Xs\h=o — m/h where h is a small applied 
magnetic field of the order J/100 and m is the corre- 
sponding magnetization which can be calculated directly 



FIG. 4: Leading correlation lengths at the S'l7(4)-point. The 
symbols at T — are CFT results where the short lines should 
indicate that these results are expected to be valid (up to 
log. corrections) also at small but finite temperatures. The 
inset shows the temperature dependence of the wave vector k 
in the case of incommensurate oscillations. 

in the asymptotic expansion (Eq. JJ5J) of the dimer corre- 
lation functions {SqSIS*S* +1 ) and (tq rf r^r^ +1 ). £d is 
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in T/J £ [0.2 : 0.5] only slightly smaller than the largest 
CL supporting the picture of a short range dimcr order 
in this temperature regime. For lower temperatures T£d 
starts to decrease as is expected because the ground state 
is not dimerized and £d should therefore remain finite. 
Note that due to the SU(4) symmetry the relation 

(S?S?) = (T?Tf) = 4<Sf rf SjVf ) (8) 

is valid for the thermodynamic correlation functions with 
arbitrary components a, (3 = x, y, z, 1 - 6 Therefore the CLs 
depicted in Fig. ^ as straight and dashed lines have 
connected non-zero matrix elements for the spin-spin, 
orbital-orbital and mixed spin-orbital correlator. At high 
temperatures these CLs show commensurate oscillations, 
i.e., the wave vector k is equal or ir. The 2 CLs which 
show 7r-oscillations at high temperatures (straight lines 
in Fig. 0J) merge at a temperature T/J ~ 0.5 with a 
square root singularity (note that the number of data 
points is relatively small in this region) and the oscil- 
lations become incommensurate and temperature depen- 
dent at lower temperatures as shown in the inset of Fig. 0] 
A similar behavior has been observed for the longitu- 
dinal spin CL in the X X Z-model for finite magnetic 
fields Note that incommensurate oscillations are associ- 
ated with a pair of complex conjugate eigenvalues of the 
transfer matrix (see Eq. J^J) which can occur because the 
transfer matrix is real but non-symmetric. Additionally 
a CL (dashed line) is present which is non-oscillating for 
all temperatures. 

It is known that the energy levels Ei of a finite chain 
of length L show a finite-size scaling behavior 



where Eq is the ground state energy and x% is the scaling 
dimension of the operator leading to the excited state 
with energy Ei $L We can get the corresponding result for 
an infinite chain at finite temperature by the substitution 
L-fv/T. This leads to 

er^lnAo-lnlA.H^ (10) 

v 

where Aj, Ao are eigenvalues of the transfer matrix. 
The possible scaling dimensions for an SU(N) symmet- 
ric model are given by x p — p(N — p)/N with p = 
1,2, ■ ■ ■ , N — 1. The dominant mode belongs to p = 1 
and oscillates at zero temperature with k = 2tt /N 
In the case of SU(4), scaling dimensions 3/4, 1,3/4 are 
therefore possible and the dominant mode has x = 3/4 
and oscillates with k = ir/2. Using Eq. I|10|) we obtain 
£„./2 = 1/12T. The non-oscillating CL must belong to 
the second possible scaling dimension x = 1 leading to 
£o = 1/16T in the low-T regime. However, this 1/T be- 
havior will be directly observable only at very low tem- 
peratures (see Ref . |l!j) . At the lowest temperatures we 
are able to access in this study it is covered by logarithmic 
corrections which are present here due to marginally ir- 
relevant operators in the low-energy effective theory^i^ 



Cardy has shown that such logarithmic corrections have 
universal character and can be displayed in leading order 

as 

2-kv ( 2b n 1 \ 

where b n , b are universal numbers^ However, it should 
be stressed that this asymptotic form is only observable 
for very large chain lengths L when the logarithmic cor- 
rection becomes independent of the coupling constant g 
of the considered marginal operator. Otherwise we can 
only replace g by its renormalization group improved 
value g(lnL) so that we must substitute InL — > \n(L/Lo) 
in Eq. I|ll|) with a non- universal Lq. But even in this case 
Eq. i|ll[l remains useful because it describes the scaling 
behavior of different excited states by a single parameter. 
By L -> v/T we find for the CLs 

tl ~ v [ l+ b bxT /Tj [U) 

where T is non-universal. This means that logarithmic 
corrections lead to an effective temperature dependent 
correction of the scaling dimension xi. The universal 
factor 2b n /b for the SU(4) model has been derived in 
Ref. l35l and is given by 2b n /b = —1/16 for scaling dimen- 
sion x = 3/4 and by 2b n /b = — 1/4 for x = 1. Therefore 
the logarithmic correction for £0 is more important than 
that for (,k/2 i n qualitative agreement with the data. To 
illustrate more quantitatively that the logarithmic cor- 
rections can explain the deviations from the expected 
CFT values, we have estimated To at T — 0.05J and 
T = 0.0625 J for £ n/2 by applying Eq. (T2J and used this 
value to calculate £0 at the same temperature. The re- 
sults are depicted as plus-symbols in Fig. ^ and are in 
good agreement with the numerically calculated We 
like to remark that this is a very sensitive test because 
To depends exponentially on £. Note also that we cannot 
extend this comparison between field theory and numer- 
ics to larger temperatures because we already know that 
the low energy effective theory breaks down at least at 
T/J w 0.08. 

IV. GAPLESS PHASE (5J7(4) WITH RESCALED 
VELOCITIES) 

By an RG analysis it has been shown that the Hamil- 
tonian J2J is not only critical at the 5'[/(4)-symmetric 
point but rather exhibits an extended region in phase 
space around this point which is gapless (phase V in 
Fig^)ii£i2£ Considering the parameters x = 1/4 + Sx, 
y = 1/4 + 5y in Eq. J2Jl, where Sx, 5y are small perturba- 
tions, the findings can be summarized as follows: (1) For 
Sx = Sy < the breaking of the SU(A) symmetry down 
to SU{2) x SU{2) is irrelevant and the low-energy be- 
havior is still described by the 5C/(4)i WZW model. (2) 
Sx + 5y < with 8x ^ Sy is marginally irrelevant and the 
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critical theory is the SU{2) 2 x SU(2) 2 WZW model with 
unequal velocities v s , v a , but with unchanged scaling di- 
mensions. (3) Sx + Sy > is marginally relevant leading 
to a Kosterlitz-Thouless type transition into a dimerized 
phase (phase IV in Fig^l. 

Representatively for the gapless phase with unequal 
velocities of the elementary excitations we will study the 
point x = 0.1, y = 0.25. Due to the different velocities 
v s 7^ v relation JHJ is no longer valid and the degeneracy 
between spin-spin, orbital-orbital and mixed spin-orbital 
CLs is lifted. This makes the situation far more complex 
than in the St/ (4) case with numerous crossovers between 
different CLs occurring in the non-universal regime at el- 
evated temperatures. To keep it concise we show only a 
selection of CLs at high temperatures in Fig. [5] and show 
the 6 largest CLs at the lowest numerically accessible 
temperatures separately in Fig. El Over a wide temper- 




T/J 



FIG. 5: Some of the leading spin and orbital correlation 
lengths for (x,y) = (0.1,0.25). The wave vectors are shown 
in the inset if incommensurate oscillations occur. The circles 
in the inset denote ±7r/2. 

ature range a 7r-oscillating spin CL is leading. This is 
expected because antiferromagnetic spin correlations are 
stabilized with decreasing x. As in the SU (4)-symmetric 
case two CLs with 7r-oscillations merge at a tempera- 
ture T/J ~ 0.4 and show incommensurate oscillations 
below which approach tt/2 at T = 0. However, in the 
SU(4) case these CLs are threefold degenerate belonging 
to spin, orbital and mixed spin-orbital correlation func- 
tions. Here no degeneracy is present and the CLs only be- 
long to the mixed spin-orbital correlator. The orbital CL 
which shows tt/2 oscillations at T = develops incom- 
mensurability already at higher temperatures as shown 
in the inset of Fig. [5] Hitherto we have ignored the dimer 
CL £d which is not only one of the largest CL as in the 
SU(4) case but is indeed leading in a certain temperature 
range. This is a little bit astonishing because we have 
moved away from the dimer phase by Sx < 0. It shows 
that the spin-orbital model has an intrinsic tendency to- 
wards dimerization at finite temperature even in regions 
of phase space where the ground state is not dimerized. 



For a system consisting of weakly coupled spin-orbital 
chains this suggests that long range dimer order can be 
stabilized at finite temperature and that such a system 
may show a phase transition from a dimer phase in this 
temperature regime to a non-dimerized phase at lower 
temperatures. In the limit T->0we expect that £d re- 
mains finite and Fig. [5] shows indeed a decreasing T£d . 
However, even at the lowest numerically accessible tem- 
perature £d is still the second largest CL. Dominating at 
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FIG. 6: Leading spin and orbital correlation lengths for 
(x,y) — (0.1,0.25) at low temperatures. 

low temperatures is a non-oscillating spin CL. Further- 
more we see a crossover between the mixed spin-orbital 
and the orbital-orbital CL with incommensurate oscilla- 
tions with the orbital-orbital CL becoming the larger one 
at low-T. From CFT we also expect a spin-spin CL with 
incommensurate oscillations which seems to be missing 
in our data. We expect that the 2 spin-spin CLs with 
k = it depicted as dashed lines in Fig. will merge at 
lower temperature yielding this missing CL. This implies 
that we still have not reached the temperature regime 
where the behavior of the CLs becomes universal and 
can be described by CFT. Therefore we cannot compare 
the numerics with CFT in detail. However, assuming 
that we are close to the conformal regime and ignoring 
logarithmic corrections completely we can extrapolate all 
CLs (except of the dimer CL) shown in Fig. starting 
from the numerical data at the lowest accessible temper- 
ature according to their expected 1/T behavior. This 
allows us to show that the data are in rough agreement 
with CFT. 

The non-oscillating spin-spin CL must belong to an 
excitation with x s = 1 so that £ = v s /2ttT leading to 
v s 0.5. For the corresponding non-oscillating orbital- 
orbital CL we find numerically (data not shown) £(T = 
0.047) re 0.86 leading to v a m 0.25. We expect that 
the ST/ (4) representation is decomposed symmetrically 
into a spin and orbital sector both with central charge 
c = 3/2 and fundamental fields with scaling dimensions 
x = 3/8. The only difference between the two sectors 
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are the velocities. All CLs with k — tt/2 at T — are 
therefore given by an excitation with x s = x a = 3/8 so 
that we expect from CFT 
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8 v v s 



J_)T 



0.071 

T 



(13) 



where we have substituted the velocities v s , v which 
we have determined from the non-oscillating spin and 
orbital CL. At finite temperatures the mixed, spin 
and orbital CL with incommensurate oscillations dif- 
fer from each other because each has different logarith- 
mic corrections »ii We notice that the value estimated in 
Eq. (|13|l is indeed in rough agreement with the numeri- 
cal data at the lowest accessible temperature (see Fig.|SJ|. 
The mixed spin-orbital CL with k — is given by an ex- 
citation with x s — x — 1/2 leading to 
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27r i(J_ + i)T 

2 v v s v s 



0.053 
T 



(14) 



predicting §p°(T = 0.047) w 1.12 which is in 
good agreement with the numerically calculated value 
£T°(T = 0.047) re 1.09 (data not shown). 

The CFT result for the free energy of the SU(4)- 
symmetric model (Eq. (0) can be easily generalized to 
the case of SU(2) x 5'[/(2)-symmetry with different ve- 
locities 



/ = e -f (i + i]7- 

6 \ D, v a 



(15) 



where now c = 3/2. In Fig. [7|the numerically calculated 
free energy and the low-T asymptotics according to this 
formula are shown. Here we have used the spin and or- 
bital velocities which we have evaluated from the CLs and 
eo as a fitting parameter. We find eg ~ — 0.2197J and a 




In the inset of Fig. the entropy and the low-T (S ~ 
3ttT) and high-T asymptotics (S ~ In 4) are shown. Here 
the low-T asymptotics is determined solely by the spin 
and orbital velocities without any additional fitting pa- 
rameter showing that CFT for the SU(2) x SU(2) model 
with our velocity estimates gives a consistent description 
of all thermodynamic quantities at low energies. Next 
we want to consider the specific heat C shown in Fig. [S] 
Contrary to the 5'[/(4)-symmetric case not only a peak 
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FIG. 8: Specific heat (solid line) and low-T asymptotics 
(dashed line) C ~ 3ttT at (x,y) = (0.1,0.25). 

and a shoulder but instead two clearly distinguishable 
peaks are visible. Note that due to the different veloci- 
ties, orbital and spin excitations have now different en- 
ergy scales. However, any saturation of one type of ele- 
mentary excitations (spinons or orbitons) will inevitably 
destroy also a description in terms of these elementary 
excitations for the other sector due to the strong cou- 
pling between both. That this is indeed the case is ob- 
vious from the susceptibility data shown in Fig. El which 
show that spin and orbital susceptibility are peaked ap- 
proximately at the same temperature^ The second peak 
of C is in a temperature range where short range dimer 
order exists and may be attributed to excitations con- 
nected to this. Finally we want to compare the suscep- 
tibility data with CFT. As in the SU (4) case we expect 
Xs,o = k/2-Kv s>0 with topological coupling constant k = 2, 
however, now v s ^ v a . Using the velocity estimates from 
the fit of the CLs we find Xs ~ 2/-7T and \ Q pa i/ir pre- 
dicting \o at T — being twice as large as \s ■ Although 
we cannot access numerically temperatures low enough 
to show that the susceptibilities converge to these values 
as T — > we note that \o is indeed more than twice as 
large as Xs at low temperatures. 



FIG. 7: Free energy and entropy (inset) at (x,y) — (0.1,0.25) 
with low and high-T asymptotics (dotted lines) as described 
in the text. 

good agreement with the numerical data for T/J < 0.1. 



V. DIMERIZED PHASE 

As described in the previous section perturbing the 
S't/(4)-symmetric Hamiltonian by Sx + Sy > is 
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FIG. 9: Spin (main figure) and orbital (inset) susceptibili- 
ties at (x,y) = (0.1,0.25). The dots denote CFT results as 
described in the text. 
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FIG. 11: Entropy and specific heat (inset) at (x,y) = 
(0.5,0.5). The dashed lines denote low-T fits as described 
in the text. 



marginally relevant. As a result a gap A in the exci- 
tation spectrum is expected. Here we do not want to 
investigate this phase transition in detail but rather re- 
strict ourselves to study the point x — 0.5, y — 0.5 as 
a representative for the dimerized phase. First, con- 
sider the spin susceptibility Xs shown in Fig. 1101 A 



* 0.3 




FIG. 10: Susceptibility (solid line) for {x,y) = (0.5,0.5) 
where the dotted line denotes the best low-T fit with S = 
0.491 exp(-0.09/T)/v / T. 

sharp peak at T = 0.12J is visible. At lower temper- 
atures \s goes to zero as expected for a gapped system. 
The usual quadratic dispersion of elementary excitations 
above the gap yields in ID an expected low-T asymp- 
totic behavior Xs ~ exp(— A/T)/y/T. Using this for- 
mula for a fit of the data in the low-T regime we find 
A = 0.090 ± 0.005 where we get a certain error estimate 
by varying the fit interval^! This value is in agreement 
with a recent series expansion study^i In Fig. II li the cor- 
responding entropy S is shown. For the free energy we 



expect at low temperatures / ~ — T 3 / 2 exp(— A/T) lead- 
ing to S ~ (AVT + B/VT) exp(-A/T). Using A, B, A 
as parameter in a fit of the numerical data for T < 0.1J 
we get the result shown as dotted line in Fig. ^] which 
is an excellent fit even up to T ~ 0.25J. From this fit 
we find A = 0.094 ± 0.004 where the error stems again 
from a variation of the interval used in the fit. We note 
that this value is in perfect agreement with the estimate 
from the susceptibility fit. In the inset of Fig.^Jthe spe- 
cific heat (solid line) is shown. The dotted line denotes 
the low temperature specific heat calculated from the en- 
tropy fit and shows good agreement with the numerical 
data. Note that the excitation energy 8 for a single bond 
is now J/4 so that the peak in the specific heat at ap- 
proximately 8/2 can be attributed again to excitations 
of the chain with dimer order which will now become a 
true long range order at T = 0. Surprisingly we still 
find a shoulder in C at higher temperatures as in the 
SU(4) case. Comparing with the CLs plotted in Fig. 1121 
we find indeed that the eigenspectrum of the transfer ma- 
trix at elevated temperatures has not changed drastically 
compared to Fig. Especially, we find again a dimer 
CL (dot-double dashed line) which becomes the second 
largest CL in T/Je [0.3 : 0.55] which coincides with the 
shoulder in C. Whereas this dimer CL stays finite as 
T — > another dimer CL (also depicted as dot-dashed 
line)22i is diverging more strongly than 1/T indicating 
long range dimer order at T = 0. Because x = y the 
CLs belonging to the spin and orbital correlation function 
are equal, however, CLs belonging to the mixed correla- 
tor are different in general^ The two CLs (solid lines) 
which merge at T / J ~ 0.5 yielding a CL with correspond- 
ing incommensurate oscillations below that temperature 
have indeed only non-zero matrix elements for (SqS*) 
and (tqT^). The 7r-oscillating CLs with non-zero matrix 
element for (S§t§S*t*) merge already at a higher tem- 
perature and the oscillations become incommensurate be- 
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FIG. 12: Leading CLs at (x,y) = (0.5,0.5). The inset shows 
the dependence of the wave vector on temperature in the case 
of incommensurate oscillations. 

low as shown in the inset of Fig. 21 Similarly, the largest 
non-oscillating CL belongs only to the spin and orbital 
correlator whereas the second largest appears only in the 
asymptotic expansion of the mixed correlator. The nu- 
merical data indicate that all CLs apart from the dimer 
CL stay finite as T — > which is expected for a system 
with a gap. 

VI. CONCLUSIONS 

In summary, we have used the TMRG method to study 
thermodynamic properties of a spin-1/2 SU(2) x SU(2) 
spin-orbital model. We have concentrated on three rep- 
resentative points in the phase diagram: (1) The SU(i) 
symmetric point where the model is Bethe-ansatz in- 
tegrable and shows critical properties at T = 0, (2) a 
point where the model still exhibits critical behavior at 
zero temperature, however, the symmetry is reduced to 
SU(2) x SU(2), and (3) a point where the ground state 
has long range dimer order. 

In the first two cases we have compared our numerical 
results in detail with CFT in the low temperature regime 
and with special attention to the evolution of various CLs 
with temperature. In both cases three different regimes 
could be identified: At elevated temperatures crossovers 
between different CLs are frequent and the correspond- 
ing wave vectors are commensurate. The 7r-oscillating 
CLs belonging to the spin, orbital or mixed spin-orbital 
correlator merge at a certain temperature with a square 
root singularity and the oscillations become incommensu- 
rate and temperature dependent below which character- 
izes the intermediate temperature regime. Finally, at low 



temperatures the conformal invariant regime is reached 
where the CLs belonging to critical excitations diverge as 
1/T up to logarithmic corrections and no crossover oc- 
cur. The comparison of the numerics with CFT was espe- 
cially simple at the ST/ (4) point because the velocities of 
the elementary excitations are known exactly from Bethe 
ansatziiS In this case we have also been able to show 
that the predicted logarithmic correctionsiZ^ account 
quantitatively for the observed deviations from a pure 
£ ~ 1/T behavior at low temperatures. In the second 
case a comparison with CFT was far more complicated 
as the velocities v s , v a are a priori unknown and several 
crossovers in the lowest numerically accessible tempera- 
ture region show that the conformal invariant regime was 
not reached. However, we have been able to get rough 
estimates for the spin and orbital velocities from the CLs 
and have shown that these estimates give a consistent de- 
scription of other thermodynamic quantities within CFT. 

For the point in the dimerized phase we have calcu- 
lated the gap from the susceptibility and entropy data 
and have found good agreement with a recent series ex- 
pansion study^i More important, we have found that the 
behavior of CLs in the high and intermediate tempera- 
ture regime here is quite similar to the gapless phase. 
At low temperatures a dimer CL diverges more strongly 
than 1 /T whereas all other CLs stay finite as is expected 
for a system with long range dimer order at T = 0. 

Most surprisingly we find in all 3 cases a strong ten- 
dency towards dimerization in a certain finite tempera- 
ture interval and dimerization could be the leading insta- 
bility in this region even if the ground state is not dimer- 
ized. This behavior is very similar to what was observed 
recently in another spin-orbital model with spin S = 1— 
indicating that such instability may be an intrinsic prop- 
erty of these kind of Kugel-Khomskii Hamiltonians. In a 
system of weakly coupled spin-orbital chains this insta- 
bility may turn into true long range dimer order at finite 
temperature. More generally speaking, it indicates that 
even a pure electronic model may show various tempera- 
ture driven phase transitions between different spin and 
orbital orderings due to the large number of nearly de- 
generated states accessible to systems with unquenched 
orbital degrees of freedom. 



Acknowledgments 

I thank A. Kliimper and J. Oitmaa for useful discus- 
sions. This work has been partially supported by the 
DFG in SP1073 and by the Australian Research Coun- 
cil. The numerical calculations were partially performed 
using the APAC computing facilities. 



* sirker@phys.unsw.edu.au 



K. I. Kugel and D. I. Khomskii, Sov. Phys. Usp. 25, 231 



9 



(1982). 

2 Y. Tokura and N. Nagaosa, Science 288, 462 (2000). 

3 Y. Tokura, ed., Colossal Magnetoresistive Oxides (Gor- 
don and Breach Science, New York 2000), and references 
therein. 

4 J. M. D. Coey, M. Viret, and S. von Molnar, Adv. Phys. 
48, 167 (1999). 

5 G. Khaliullin and S. Maekawa, Phys. Rev. Lett. 85, 3950 
(2000). 

6 B. Keimer, D. Casa, A. Ivanov, J. W. Lynn, M. v. Zim- 
mermann, J. P. Hill, D. Gibbs, Y. Taguchi, and Y. Tokura, 
Phys. Rev. Lett. 85, 3946 (2000). 

7 J. Hemberger, H.-A. K. von Nidda, V. Fritsch, J. Deisen- 
hofer, S. Lobina, T. Rudolf, P. Lunkenheimer, F. Licht- 
enberg, A. Loidl, D. Bruns, et al., Phys. Rev. Lett. 91, 
066403 (2003). 

8 M. Cwik, T. Lorenz, J. Baier, R. Miiller, G. Andre, 
F. Bouree, F. Lichtenberg, A. Freimuth, R. Schmitz, 
E. Miiller-Hartmann, et al., Phys. Rev. B 68, 060401 
(2003). 

9 M. Mochizuki and M. Imada, Phys. Rev. Lett. 91, 167203 
(2003). 

10 G. Khaliullin, P. Horsch, and A. M. Oles, Phys. Rev. Lett. 
86, 3879 (2001). 

11 Y. Ren, T. T. M. Palstra, D. I. Khomskii, A. A. Nugroho, 
A. A. Menovsky, and G. A. Sawatzky, Phys. Rev. B 62, 
6577 (2000). 

12 Y. Ren, T. T. M. Palstra, D. I. Khomskii, E. Pellegrin, 
A. A. Nugroho, A. A. Menovsky, and G. A. Sawatzky, Na- 
ture 396, 441 (1998). 

13 M. Noguchi, A. Nakazawa, S. Oka, T. Arima, Y. Wak- 
abayashi, H. Nakao, and Y. Murakami, Phys. Rev. B 62, 
R9271 (2000). 

14 J. Sirker and G. Khaliullin, Phys. Rev. B 67, 100408(R) 
(2003). 

15 C. Ulrich, G. Khaliullin, J. Sirker, M. Reehuis, M. Ohl, 
S. Miyasaka, Y. Tokura, and B. Keimer, Phys. Rev. Lett. 
(2003), in print. 

16 Y. Q. Li, M. Ma, D. N. Shi, and F. C. Zhang, Phys. Rev. 
Lett. 81, 3527 (1998). 

17 C. Itoi, S. Qin, and I. Affleck, Phys. Rev. B 61, 6747 
(2000). 

18 B. Sutherland, Phys. Rev. B 12, 3795 (1975). 



19 B. Frischmuth, F. Mila, and M. Troyer, Phys. Rev. Lett. 
82, 835 (1999). 

20 Y. Yamashita, N. Shibata, and K. Ueda, Phys. Rev. B 58, 
9114 (1998). 

21 S. K. Pati and R. R. P. Singh, Phys. Rev. B 61, 5868 
(2000). 

22 P. Azaria, E. Boulat, and P. Lecheminant, Phys. Rev. B 
61, 12112 (2000). 

23 S. K. Pati, R. R. P. Singh, and D. I. Khomskii, Phys. Rev. 
Lett. 81, 5406 (1998). 

24 H. F. Trotter, Proc. Amer. Math. Soc. 10, 545 (1959). 

25 M. Suzuki, Phys. Rev. B 31, 2957 (1985). 

26 J. Sirker and A. Kliimper, Europhys. Lett. 60, 262 (2002). 

27 J. Sirker and A. Kliimper, Phys. Rev. B 66, 245102 (2002). 

28 I.Peschel, X. Wang, M. Kaulke, and K. Hallberg, 
eds., Density-Matrix Renormalization, Lecture Notes in 
Physics, vol. 528 (Springer, Berlin, 1999), and references 
therein. 

29 X. Wang and T. Xiang, Phys. Rev. B 56, 5061 (1997). 

30 N. Shibata, J. Phys. Soc. Jpn. 66, 2221 (1997). 

31 I. Affleck, Nucl. Phys. B 265, 409 (1986). 

32 I. Affleck, Phys. Rev. Lett. 56, 2763 (1986). 

33 A. Kliimper, J. R. R. Martinez, C. Scheeren, and M. Shi- 
roishi, J. Stat. Phys. 102, 937 (2001). 

34 J. L. Cardy, J. Phys. A 17, L385 (1984). 

35 K. Majumdar and M. Mukherjee, J. Phys. A 35, L543 
(2002). 

36 J. L. Cardy, J. Phys. A 19, L1093 (1986). 

37 W. Zheng and J. Oitmaa, Phys. Rev. B 64, 014410 (2001). 

38 Within this paper we will use e = 0.05 in all calculations. 

39 Note that spin and orbital susceptibility are equal along 
the line x = y. 

40 More precisely, \o is peaked at T ~ 0.043J and \s at 

0.033J. 

41 This error does not include systematic errors due to the 
truncation of the Hilbert space which are difficult to esti- 
mate. 

42 We use the same symbol for both CLs because they are 
of the same type. However, they are connected to different 
eigenvalues which are crossing at a certain temperature. 

43 Eq. © is not valid here. 



